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I.  INTRODUCTION 

The  digital  number  representing  the  Short  Term  Averager  (STA)  out- 
put of  a DIMUS  system  must  be  requantized  to  fewer  bits  in  order  to 
limit  the  size  of  further  averagers.  It  is  desirable,  however,  to 
maintain  linearity  for  small  signal  to  noise  ratios  (reference  UJ;. 

A reasonable  compromise  between  minimizing  bit  requirements  and  main- 
taining linearity  is  to  requantize  the  STA  data  so  that  it  is  approx- 
££*  linear  in  tb.  ranea  IV.  (»her.  Vk  is  the  fAtutputatandard 
deviation  for  noise  only  at  theTiydrophones)  and  hard-limited  for 
numbers  outside  this  range.  Therefore,  it  is  necessary  for  design 
nproosos  to  calculate  % . This  number  is  usually  found  under  the 
assumption  that  the  time  sampled  inputs  to  the  STA  are  uncorrelated. 

In  this  report  the  effect  of  the  correlation  of  the  STA  innuts  on  Tj* 
is  determined.  Computer  simulation  results  are  included  to  check  the 

calculations.  aa  v ..  . jl  / 4 


T. r 

Figure  one  is  a block  diagram  of  the  system  being  considered.  The 
probability  distribution,  mean,  variance  and  range  of  the  data  at  each 
point  have  been  derived  by  Penn  and  Trudell  (reference  (2))  and  are 
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included  here  for  convenience.  These  numbers  are  derived  as  follows s 

Assuming  the  x^t)  are  uncorrelated,  z has  the  binomial  distribution 
(reference  (3)). 

Pz(4)  = 

where  p = PrCx^o)  = l/2  for  noise  only. 

Then  _ 

2«  MP-^M  - O 

v*/=  mp(i-p)=  ~ 

Now,  assuming  M is  sufficiently  large  so  that  pz(Z)  may  be  considered 
Gaussian 


an(j  see  appendix  A. 

r TT-  Z rrl  - IT-  Z.  M 
s TT  * * IT  1 . 

If  the  space  averager  implements  the  equation 


^ * *<8;-  42 

we  see  that 

V*  = ItTf-  UBiEBi  *1  I 6£; 

^ i i S 

and  if  the  beams  are  uncorrelated 

Vi*  - I6<7w  - zoV'J 
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Finally,  if  the  are  uncorrelated 


K * K C = SK 


IT-  t 


Where  K is  the  number  of  samples. 

If  (-1,  +1)  clippers  are  used,  rather  than  (0,  1)  clippers  as  we  have 
assumed  equation  (l)  should  be  multiplied  by  4 to  find  7^, 

The  analysis  has  assumed: 

1.  Noise  only  at  the  hydrophones. 

2.  The  noises  at  the  hydrophones  are  uncorrelated. 

3.  M is  sufficiently  large  so  that  pz  (Z)  is  Gaussian. 

4.  The  beams  are  uncorrelated. 

5.  The  are  uncorrelated. 

Assumptions  1 and  3 are  reasonable;  assumption  2 is  unrealistic  for 
small  arrays  and  low  frequency  bands;  assumption  4 is  not  true  since 
the  beams  are  correlated  to  some  degree.  Assumption  5 is  the  problem 
we  now  consider.  Penn  and  Trudell  also  consider  the  problem  of  corre- 
lated samples,  but  from  a different  viewpoint. 


III.  EFFECT  OF  CORRELATION 


b = <*i 


rkl  = Tl 

l K K _____ 

s.2  I <*.i 

• 4*  I 4 - 1 
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but 

1 -T-  * 

and 

Thus 

SI3.J  - 

% 

cx.ja.i  *ince 

. , *-1 

Pa(t)  is  even. 

K 

rkl=  2 zff.lt i-iuTj 

• K«C*  * (z) 

l*  l 


where  (*a(T)  is  the  normalized  autocorrelation  function  of  a,  and  AT 
IS  the  sample  interval. 


Thus,  the  problem  is  reduced  to  finding  Pa(t)  in  terms  of  P x (t) , the 
input  autocorrelation. 

Assuming  the  beams  are  uncorrelated 

Pa^  = Pw(t) » 

so  we  must  find  pw(t)  in  terms  ofpz(t). 


Papoulis  shows  (reference  (4) , p.  482) , that  if  z is  a Gaussian 
random  variable  with  0 mean,  variance  7"jl  and  autocorrelation 
function  Pz(t)  that 


✓ 


so 


R(f)  a {tl-P.V)]*  ♦ f.tt) 

formalizing  Ry(t)ve  obtain 

V A ts  T,‘ 

= ^ fci  - f.'ftfls  *[p,(t)]  - lj  _ 


Finally  vJt  *»TC 


Ci  it)  * £ sin 

if  the  x^are  Gaussian  variables  witli  iJeet»<*l 
Thus  our  final  result  is 

e.tt)  i •»'«><(♦)* 

* LJi"  Vf  ii«’V»ct))J  - ij 
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Equation  2 hes  been  derived  under  assumptions  1,  2,  3 and  4 of 
Section  II,  and  in  addition  we  now  require  the  x*(t)  to  be  samples 
from  a Gaussian  process  wiTH  a»T*cM*el*TiM  fwearieaa. 

IV.  EVALUATION  FOR  A SPECIFIC  CASE 

It  now  remains  to  evaluate  the  effect  of  all  this  for  a specific 
case  of  interest.  Let  the  power  spectral  density  of  x,  S^f) , be  given 

Cf ) * -^5-  fe  ’ 7 < I * I < * f* 

IB  * • 


= 0 elsewhere. 

then 

o (t)  s s,n  (ire  f)  coi  ( Mrf.t) 

* ire  t 

Appendix  B lists  the  values  jf  f*(AT)  for  AT  = 40  ysec,  fQ  = 2500  Hs., 
B = 3000  Hz.  Mow  if  M = 48,  K = 512  we  obtain,  for  the  uncorrelated 

v,  s Js  (M8)(SU)  k = Zll  l«e«t  |. 

V b Y J.H 

Then  using  the  computed  values  of  the  correlation  from  appendix  Bj 
equation  2 ex.  a A J V*  fiad 

zn[i 
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Thus  the  false  assumption  that  the  a^  are  uncorrelated  causes  an  error 
of  roughly  5056  in  the  computed  value  of  the  STA  output  standard  devia- 
tion. 

V.  COMPUTER  SIMULATION  RESULTS 

R.  L.  Gordon’s  DIMUS  simulation  program  (reference  (5))  was  used 
to  check  the  calculations.  The  program  was  run  using  uncorrelated, 
Gaussian  noise  at  the  inputs  having  a +3  db/  octave  slope  in  the  1-4 
band.  Ho  space  averaging  was  used  in  order  to  eliminate  the  effects 
of  beam  to  beam  correlation.  Thus,  the  program  conforms  to  assumption 
1,  2,  3 and  4 of  Section  II. 

The  program  used  (-1,  +l)  clippers  and  K = 128. 

Using  these  parameters,  equation  1 predicts 

=•  [ (128)  ^ H7.Z2 

if  the  STA  inputs  are  uncorrelated.  From  equation  2 an  i eqvanffl  3 

rk  = *±, < «•?*)]**  w.7 s 

the  result  obtained  when  the  correlation  is  accounted  for.  This 
result  was  computed  using  the  first  20  values  of  px(AT)  for  the 
+3  db/  octave  noise  used  in  the  program. 

Table  one  shows  the  standard  deviation  of  the  STA  output  for  all 
48  beams  as  computed  from  the  simulation  program.  These  standard 
deviations  were  each  computed  from  129  data  points.  The  average  of 
this  data  is  75.57  and  the  data  has  a standard  deviation  of  7.8.  This 
number  is  close  to  the  result  obtained  from  equation  2,  and  we  conclude 
that  the  computer  simulation  verifies  the  result  of  the  analysis. 

VI.  CONCLUSIONS 

^It  has  been  shown  that  the  numerical  value  to  use  for  the  standard 
deviation  of  the  STA  output  is  roughly  1.4  times  the  number  calculated  - 
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assuming  uncorrelated  samples,  if  assumptions  1-4  are  true.  Computer 
simulation  results  verify  the  analysis.  The  effect  of  beam  to  beam 
correlation  and  of  correlated  noise  at  the  hydrophones  on  the  STA 
output  standard  deviation  has  not  yet  been  determined,  and  is  the 
subject  of  further  study,  ^ 
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APPENDIX  A 
Derivation  of  u}  and 

Assume  is  a Gaussian  random  variable  with  zero  mean  and 
variance  , then  if 

w - |i| 

_ rW 

[w  i Wj  s 2.  Pp*b  ^2  £ Wj  * £ J 

o 

then  K (w)  = L p„t  [ W £ Wj  = Z )**  l W) 

t - —l 

=1fIFvi  e ZVt  wi0- 

A /ou>  u/e  *f irt J w » 

u5  - / ^ pw(w)Ju>  * (n?  ~ e 

— «0  o"  1 

let 


d t - u/  J u)  J wJw  1 ^ d i 
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Next,  we 

•f»ad  Vwlj**>kerte 

1 * 

II 

tv*  - 

- Cu;)2. 

Vow 

let 

f„7F 

u/ 

e &M*i 

t = ^ 

frvi 

1*ke 

n 

and 

[ tv1  - IVi*  t 

so 


ool  Ju>  * if?  V.Vjt 


and 


1 * flJt 

YIP  Vo 


(A/ 


VtP(I) 

but  the  Integral  is  equal  to 

— > 1 

<*'  = v*\ 


so 

Therefore 


Vwl=  u.*-(£DJ  = v," 
V*\ 


_ TT-2.  — i 


TT 


‘Reference  Data  for  Radio  Engineers,  ITT,  4th  edition,  p.  1079, 


egv.  144. 
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PCf)  * - cos(iff.t) 

7TB  t 

UHOyStc)  = ^.<p(TTt3<ll>1J(.»<M»l)~lJ  (<>UT(tii)t(H|3i) 

Tf  (3Xioi)Coi  x\o~J) 

» 

1 

= .in 

PCTOi  * .It  | 

P(Ui>)=  -.1*0 

- 

P(U«*  -.534 

p(2ooj  = -.505" 

P(1H0)  = Z7  C 
PtitO)  = .osi 
P(32  (?)  = -Oil 
PCJtOj  i-.Ofc 
PtHftO)  s.  ,|U 
PC  •/Ho)  = IH1 

e(n*o)  = -.043 

PiSlO)  - -.043 

K I 
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